Relativistic Hartree-Fock theory. Part I: density-dependent effective Lagrangians 
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. Effective Lagrangians suitable for a relativistic Hartree-Fock description of nuclear sys- 

O ' 

£SJ ' terns are presented. They include the 4 effective mesons o~,ui, p and ir with density-dependent 



meson- nucleon couplings. The criteria for determining the model parameters are the repro- 



duction of the binding energies in a number of selected nuclei, and the bulk properties of 
nuclear matter (saturation point, compression modulus, symmetry energy). An excellent 
description of nuclear binding energies and radii is achieved for a range of nuclei encompass- 
^ ' ing light and heavy systems. The predictions of the present approach compare favorably 

with those of existing relativistic mean field models, with the advantage of incorporating the 

o ' 
oo 
O 

)Q ■ PACS numbers: 21.30.Fe, 21.60.Jz, 21.10.Dr, 24.10.Cn, 24.10.Jv 

O 



effects of pion-nuclcon coupling. 



£ ' I. INTRODUCTION 



The most actively investigated fields in present days nuclear structure and reaction studies are the 
$— i ' production and exploration of iso 

niaaaaa 



exotic nuclei 



oratio_n of isotopes with extreme neutron-to-proton number ratios, the so-called 
. In the last 20 years, technological breakthroughs in producing 
opened up entirely new and exciting frontiers for nuclear physics P]. 
With the development of unstable nuclear beams 

p,a 

,a number of unexpected phenomena have 
seen discovered like neutron and proton skins and halos Ifll lit 



iL u s 

3, 



been discovered liKe neutron and proton sKms and naios in. im. nil. 11211. modifications of shell closures 

matter distributions |16l. 1171] . etc. With further developments, many other new features are likely to 
be found. 

In the universe, novae and supernova explosions leading to neutron stars, X-ray and gamma-ray 
bursts, all depend upon reactions involving nuclei that do not naturally occur on Earth. There 
still remain many mysteries in the origin of elements with the rapid neutron-capture process 

accounting for the formation of about half of A > 60 stable nuclei in nature. The r-process is the most 
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complex nucleosynthesis process from the astrophysics as well as nuclear physics point of view |19L |2C ] . 
Understanding the time scales and energies for such processes, and hence the stellar evolution itself, 
should be connected with the understanding of the exotic nuclear systems, either in astrophysical 
conditions or in the extreme nuclear environment. The exotic nuclei thus provide also a new testing 
standard for our understanding of phenomena of astrophysical interest. 

Many of the above physics issues can be successfully studied in the framework of self-consistent 
mean field approaches. Non-relativistic Hartree-Fock approaches based on Skyrme type forces I21II22I . 
local density approximations of the Brueckner G-matrix 23, 2^ |. and effective Gogny type forces 



have been developed. A quantitative understanding of nuclear matter as well as finite nuclei was thus 
obtained. In parallel, the relativistic Brueckner-Hartree-Fock approach was applied to nuclear matter, 
at first without full self-consistency [3] , and then fully self-consistent calculations were carried out by 



Brockmann and Machleidt 



27 



;t 29]. I 

interaction and combined with relativistic Brueckner-Hartree-Fock approach, the relativistic approach 
offers an appealing framework to investigate and explore macroscopic properties of hot and dense 
nuclear matter |29(. Using a density-dependent parametrization of the Dirac-Brueckner G-matrix 
in nuclear matter, an effective Dirac-Brueckner-Hartree-Fock model for finite nuclei was achieved by 
Boersma and Malfliet which also offers a wide range of applications in the nuclear physics 

and astrophysics domains. 

During the past years, the relativistic mean field (RMF) theory 32, 331 has received wide attention 

nnnn 

due to its successful description of numerous nuclear phenomena |34L 13.1 I361 . 13 71 13SL 1391 ] . The RMF 

can be viewed as a relativistic Hartree approach with a no-sea approximation. In the framework of the 
RMF, the nucleons interact via the exchanges of effective mesons, and photons. For nuclear structure 
applications, the relevance of relativity is not the need for relativistic kinematics but that of a covariant 
formulation which maintains the distinction between scalar and vector fields (more precisely, the zeroth 
component of Lorentz four- vectors) . The representations with large scalar and vector fields in nuclei, 
of the order of a few hundred MeV, provide simpler and more efficient descriptions than non-relativistic 



approaches 



40 



41 



421 ] . The dominant evidence is the spin-orbit splitting. Other evidence includes the 



density dependence of the optical potential, the observation of approximate pseudo-spin symmetry, 
etc. With a very limited number of parameters adjusted to reproduce selected empirical observables, 
the RMF theory can describe very satisfactorily general properties such as the masses and radii of 
nuclei in the whole periodic table, or specific aspects like the observed kink in the isotopic shifts of 



the 



44 



b-region |43j. It gives naturally the s pin- orbit potential, the origin of the pseudo-spin symm etry 



451 ] as a relativistic symmetry 
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4a ] . the spin symmetry in the anti-nucleon spectrum 



49|, 



etc. It also performs quite well for nuclei far away from the line of /3-stability with proper treatment 



of the pairing correlation and continuum effects 



3 

In spite of the success, there are still some shortcomings in the RMF theory. A too small value of 
the symmetry energy would be obtained if the coupling constant g p takes on its experimental value. 
Actually, important contributions due to the exchange terms are missing in the model. Among the 
effective mesons which mediate the nucleon interaction, the pion meson might be the most important 
one. Because of the limit of the approximation itself, the one pion exchange cannot contribute within 
the framework of the RMF if parity is to be conserved. The RMF cannot take into account the effective 
spin-spin interaction coming from the nucleon density through the anti-symmetrization of the wave 
function. This effect is more important in spin-unsaturated systems 53]. 



Early attempts to investigate the structure of nuclear matter and finite nuclei in the relativistic 
Hartree-Fock (RHF) approximation have been made The effective Lagrangians did not include 

density dependence or meson self-couplings. They gave a satisfactory description of the charge densities 
of magic nuclei but failed on the incompressibility. The nuclei were not bound enough, this defect being 
more important for light nuclei. By introducing additional density-dependent terms in the Lagrangian 
density, better ag reement with the experiment and significant improvement on the incompressibility 



were obtained 



5 a] . The density dependence was simulated by the nonlinear self-couplings of the cr-field, 



but chiral symmetry is not conserved. To recover the chiral symmetry, the non-linear self-interaction 
(NLSI) of the scalar field with zero-range was introduced in Ref. p7j |. With the zero-range limit, 
the nonlinear self-couplings are expressed in terms of the products of six and eight nucleon spinors, 
where the exchange contributions can be evaluated by the Fierz transformation. It is shown that, in 
the RHF formalism the explicitly isovector-independent NLSI generates strongly density-dependent 
scalar, vector and tensor nucleon self-energies. The description of nuclei thus obtained is in reasonable 
agreement with experimental data. 

Although some progress have been obtained in the relativistic Hartree-Fock description of nuclear 
systems, it is still a long way compared to the RMF, either on the quantitative descriptions of nuclear 
systems or their extrapolation to the exotic region. Especially compared with the wide applications 
of RMF, the disparity seems more remarkable. Of course, the technical and numerical difficulties 
are increased strongly with the inclusion of the exchange terms. On the other hand, there is still no 
appropriate effective Lagrangian for the RHF approach, which does not break the chiral symmetry 
nor increase the complexity of the theory itself. 

Based on the above considerations, we have constructed a density-dependent relativistic Hartree- 
Fock (DDRHF) theory by quantizing a Hamiltonian which contains density-dependent meson-nucleon 
couplings [s2|. For the first time, the RHF approach is systematically applied to effective Lagrangians 
with density-dependent meson-nucleon couplings. The main finding is that, with a number of ad- 
justable parameters comparable to that of RMF Lagrangians, one can obtain a good quantitative 
description of nuclear systems without dropping the Fock terms. The pseudo-spin symmetry in finite 



nuclei within our DDRHF model has been investigated in Ref. [5^. In this paper we present the 
method for obtaining the DDRHF effective Lagrangians and the results for the selected nuclei. Sys- 
tematical applications of the DDRHF in nuclear matter and finite nuclei will be shown in a forthcoming 
paper Q. 

This paper is organized as follows. The derivation of the DDRHF is briefly introduced in Section 
ITT1 In Section II I II are given the parametrizations of DDRHF and some initial applications in nuclear 
matter and in the set of selected nuclei. A short summary is given in Section IIVI The calculations of 
center-of-mass corrections for energies, radii and charge densities are explained in the Appendix. 

II. THEORY FRAMEWORK 
A. General Formalism 

The general framework can be found in several papers published in the literature. We follow closely 



the formalism and notations of Refs. 



54 



55f | . The detailed expressions for nuclear matter and finite 
nuclei in the spherical symmetry approximation were already given in these two references, except that 
they were established for the case of density-independent coupling constants. Therefore, we recall here 
only the main steps and we briefly explain how to treat the additional rearrangement terms brought 
about by our density-dependent Lagrangian. 

The starting point is an effective Lagrangian density C{x) constructed from the degrees of freedom 
associated with the nucleon (■0), two isoscalar mesons (a and lu), two isovector mesons (ir and p) 
and the photon (^4) fields. The stationarity condition of the action integral J d 4 xC(x) variations of 
the physical fields (f> (0 = tp, <r, uj, p, ir and A) leads to the Euler-Lagrange equations, from which one 
can deduce the equations of motion for the meson, photon and nucleon fields jo^) ]. The meson and 
photon fields obey inhomogeneous Klein-Gordon equations and a Proca equation with source terms, 
respectively, whereas the nucleon field obeys a Dirac equation. 

The Hamiltonian operator, i.e., the (00) component of the energy-momentum tensor T^ LV = 



d u (f)i — g^ u Jf, can be formally obtained through the general Legendre transformation 

The next step is to use the Klein-Gordon and Proca equations to formally express the meson and 
photon fields in terms of nucleon fields ip only. Thus, one would be able to obtain the Hamiltonian in 
nucleon space in a form suitable for studies of nuclear systems. We make the simplifying assumption 
of neglecting the time component of the four-momenta carried by the mesons, i.e., the meson fields 
are assumed to be time independent. This amounts to neglect the retardation effects. The energies 
involved are small compared to the masses of the exchanged mesons, so that this approximation should 



be valid for the a-, u- and p-induced interactions, and also, to a lesser extent for the pion. Then, the 
meson propagators have the usual Yukawa form, 

]_ e -m a \Xi— X2\ 



D a {x 1 ,x 2 ) 



Att \x\ — a?2 1 



(2) 



and the nucleonic Hamiltonian can be written as 



H = j d 3 x $ [-z 7 • V + M] t/,] + d 3 Xl d 4 x 2 ^ $(xi)i>(x 2 )T i (l,2)D i (x 1 ,x 2 )' l p(x 2 )iP(x 1 ), (3) 



t— ct,uj, 
P,tt,A 



where the interaction vertices are generalizations of those of Ref. 

r CT (l,2) =-g a (l)g a (2), 

r w (l,2)= + «7 w (l) 7M (l)^(2)y(2), 

T p (l, 2) = + ^(1)7^(1)^(1) • ^(2)y*(2)r(2), 



rir(l,2)=-f^7*yfe7^ 



r A (l,2)EE + T [ 7At (l-T3)] 1 [7 /i (l-r 3 )] 2 . 



The nucleon field operators ip and ijj' are expanded as 



i 

^0*0 =E [fh^M+ak^e-^d, 



(4a) 
(4b) 
(4c) 
(4d) 

(4e) 

(5a) 
(5b) 



where fi(x) and gi(x) are complete sets of Dirac spinors, Cj and cj represent annihilation and creation 
operators for nucleons in a state i, while di and d\ are the corresponding ones for antinucleons. Since we 
will restrict ourselves to the Hartree-Fock approximation level and study the exchange corrections to 
the mean field approach, we keep the same level of approximation, the so-called no-sea approximation, 
i.e., the d and $ terms are omitted in the expansions. Their inclusion leads to divergences and requires 
a cumbersome renormalization procedure |34| . The {fi} states will be determined by the self-consistent 
HF scheme. 

The trial ground state |<3?o) is the HF Slater determinant built on the lowest {/,} states. The energy 
functional, i.e., the expectation value of the Hamiltonian can be obtained as 



E = ($ | H |$ ) = <$o| T |$ > + Yl (*o| Vi |$ ) • 



(6) 



For the two-body interaction parts ($o| Vi \&o)i ° ne obtains two types of contributions, the direct 
(Hartree) and exchange (Fock) terms. For example, the u-meson contributions are: 



Ej? = - - J dx g a ($ | Ps |$o) J dx'g a ($ | Ps |$o) D a (x, x') 
E„ = + \ j dx J dx'Y [g<rUfi3] x Da(x,x') [g fpf a ] x , , 



(7a) 
(7b) 



where p s is the scalar density. The states {f a } are solutions of a Dirac equation containing a self- 
energy term E. The stationarity of the energy @ with respect to variations of the {/«} determines 
self-consistently E. 

To be complete, we would like to mention about the contributions from the one-pion exchange. 
At the Hartree-Fock level, it is well known from general arguments (chiral symmetry) and from the 
construction of the NN potential (with pair suppression mechanism) that one should adopt the pseudo- 
vector coupling to obtain reasonable results in the one-pion exchange approximation. Here, we use 
the pseudo- vector coupling for the pion, as it is done in Ref. At the Hartree (RMF) level, the 

contributions from the one-pion exchange are reduced to zero if the variational space respects parity 
conservation while in RHF these contributions appear naturally. 



In the recent years, the RMF wit 



actively developed 



61 



62 



3 



64 



i density-dependent meson-nucleon couplings (DDRMF) has been 



65 



66]. In this work, the density- dependent rclativistic Hartree- Fock 
(DDRHF) theory is derived by jointly using the RHF approach with an effective Lagrangian whose 
meson-nucleon couplings are density-dependent. Similarly to the DDRMF, the coupling constants g a , 
9uji 9p an d j v are taken as functions of the baryonic density pi,, the zeroth component of the nucleon 
current f = ip^ip. For the a- and cj-mesons, the density-dependent behaviors of g a and g u are chosen 
as 

9i(Pb) = 9i(Po)fi(0 fori = cr,w (8) 

where 

l + bi{i + d t ) 2 
1 + Cj(£ + diY 

is a function of £ = Pb/po, and pq denotes the saturation density of symmetric nuclear matter. For 
the functions fi(x), five constraint conditions = 1,/"(1) = /"(l) and /f(0) = are introduced 
to reduce the number of free parameters. 

For g p and /„., an exponential density-dependence is adopted 

9 P {Pb) = 9 P (0)e- a ^, (10a) 
fM = U{$)e- aA . (10b) 

Because of the energy-momentum conservation condition, the density-dependence in meson-nucleon 
couplings will lead to a new term in the HF self-energy, the rearrangement term E^. Thus, the total 
HF self-energy is a sum of a Hartree (direct) term Tijj, Fock (exchange) term E^ and rearrangement 
term E#. 
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B. Nuclear Matter 



In nuclear matter, the self-energy S can be written quite generally as 



(11) 



where p is the unit vector along p and pp is the Fermi momentum. The tensor piece 707 -pSy(p) does 
not appear in the HF approximation for nuclear matter. The scalar component Sg, time component 
So and space component Sy of the vector potential are functions of p = (E(p),p). 

With t he g eneral form of the self-energy, the Dirac equation in infinite nuclear medium can be 
written as I5J1 



[7 • p* + Mg] u(p, s) = 7oE*u(p, s). 



(12) 



The starred quantities are defined by 



p* =p + pY, v (p), 

M* S =M + Z S (P), 
E* =E(p) - E (p), 



(13) 



where Mg is the scalar nucleon mass, and E is the single-particle energy in the medium. The calcula- 
tion of the self-consistent self-energies follows the method indicated in Refs. 



54 



. In DDRHF, the 



density-dependence of the coupling constants gives additional contributions to the vector self-energy 
Sq, namely the direct (S^) and exchange (S^) rearrangement terms S#: 



'R 



^R,cr + ^R,lu + ^R,p + + + ^R,p + ^ 



E 

R,TT 



= s£ + s E 



R ' 



(14) 



where 



yD _ i 
^R - — 
PO 



^R — — 



m^ s dx m? i f b dx 



5e% dg a 6e t 



+ 

.E 



mf b Pb Ox 



+ 



(15a) 
(15b) 



po 5g a dx 8gu dx bg p dx 5f w dx 

with x = Pb/ Po- The rearrangement term S# can also be explicitly written in terms of self-energies. 
For example, the contribution from the cr-meson reads as 



dga_ J_ J_ 
dpb g<r 4^ 71-2 



m(p)K,s(p) + Ko(p) + Hp)K,v(p) 



p 2 dp. 



(16) 



From the self-energies, the starred functions can be determined. Finally, a self-consistent 
procedure is achieved for a given baryonic density p\, and relative neutron excess (3 = (N — Z)/A as, 




New iteration. 



(17) 



C. Spherical Nuclei 

In spherical nuclei, a single-particle state with energy E a is specified by the set of quantum numbers 
a = (a,m a ) and a = (T a ,n a ,l a , j a ), where r a = 1 for neutrons and —1 for protons. The Dirac spinor 
f a (f) in the expansion © can be written as, 



fair) 



1 / iG a (r) 
F a (r)a .r 



(18) 



where x I ( r a) is the isospinor, G a and F a correspond to the radial parts of upper and lower components, 



respectively, ^ a a ma is the spinor spherical harmonics. The Dirac spinor (|18fl is normalized according 
to 



/ drfl(r)f a (r) = f dr [G 2 a (r) + F 2 a (r)] = 1. 



(19) 



With the help of the Dirac spinors ()18|) the expression for the total energy (jEJ) can be obtained 55]. 
Then, the Hartree-Fock equations can be derived by requiring the total binding energy Eb 



Eb 



H \(f> ) - AM = E - AM 



(20) 



to be stationary with respect to norm-conserving variations of the spinors f a (or, of G a and F a ): 







5 E B -Y, E « [ flfadr 

The Hartree-Fock equations in coordinate space take the form: 

E a G a (r) = - [A _ Fa {r) + [M + S s (r) + E (r)] G a (r) + Y a (r), 
dr r 

E a F a (r) = + \ 4- + — 1 G a (r) - [M + E s (r) - E (r)] F a (r) + X a (r), 



(21) 



(22a) 
(22b) 



where Eg and So contain the contributions of all mesons to the direct (Hartree) potential and the 
rearrangement potential X a and Y a are the convolution integrals of the non-local exchange (Fock) 
potentials with F a an d G„ , respectively. The detailed expressions of the Hartree and Fock potentials 
can be found in Ref. |55j. Similarly to nuclear matter, the rearrangement potential due to the 



density-dependence of the meson-nucleon couplings can be obtained from the variation of the energy 
functional with respect to the baryonic density pb- As an example, the rearrangement potential due 
to the cr-N coupling can be written as 



4 CT) (o 



dpb 



p s (r)a(r) + 



1 1 

9a r 2 



Y,ja(Ga(r)Y^(r) + F a (r)XtHr) 



(23) 



Notice that the rearrangement potential is local. 

The system of coupled integro-differential equations 1)22(1 is rather cumbersome to solve. It is 



possible to transform it into an equivalent system of coupled differentia^ 
solved iteratively by the standard techniques used in the RMF approach 



follow the method of Ref. 



equations, which can be 



511 ] . To this end, one can 



551 ] and rewrite the functions X a {r) and Y a (r) in the form: 
F a (r)X a (r) 



_ G a (r)X a (r) _ 
Aa[r) Gl(r) + F 2 (r) Ua[r) + G£(r) + F|(r 



-F a (r) = X a , Ga (r)G a (r)+X a , Fa (r)F a (r), (24a) 



Y«(T) = ^r) ) +F2{r) Ga{r) + C^+l^l/ ^ = ^ + (24b) 
Then, the integro-differential equations (|22|) are formally transformed into a homogeneous differential 
system, 

d K a 



E a G a (r) = - 
E a F a (r) = + 



dr r 

d n a 
dr r 



F a (r) + [M + E 5 (r) + S (r) + Y ajGa (r)} G a (r), (25a) 
G a (r) - [M + E s (r) - £ (r) - ^.(r)] ^(r), (25b) 



which has the same structure as the RMF equations. This can be solved iteratively as in the RMF 
approach. 

It is very important to treat the T = 1 pairing correlations if one aims at a quantitative descrip- 



tion of finite nuclei. From a mean field point of view, one should ultimately 



Bogoliubov approach. At this stage we only use the BCS approximation 



67 



adopt a Hartree-Fock- 



681 ] to include pairing 



effects, as it is widely done in non-relativistic and relativistic mean field approaches. The pairing 



matrix elements are calculated with a zero-range, density-dependent interaction of the type 

Pb{r) 



69] 



V(n,r 2 ) = Vb^n - r 2 ) 



1 



(26) 



with Vo = -850 MeV.fm 3 and po denotes the saturation density of nuclear matter. The active pairing 
space is limited to the single-particle states below +5 MeV. The unbound states are calculated by 
imposing a box boundary condition, and one checks that the overall results are not affected by the 
choice of the box size. The Hartree-Fock-BCS equations are solved self-consistently by iteration just 



like in RMF-BCS calculations 



351 ] . This is a standard procedure where the presence of the Fock terms 



does not bring any additional difficulty. 



10 

The correction of the center-of-mass motion is treated in a microscopic way as in Ref. The 
detailed expressions are given in Appendix [A| For the radii, both direct and exchange corrections (see 
(|A19|0 are taken into account to be consistent with the theory itself. To calculate the charge densities 
and radii we take into account both the proton form factor and the center-of-mass corrections (see 
Appendix lA")). 



III. EFFECTIVE INTERACTIONS 



The long standing problem in the RHF approach is to provide an appropriate quantitative de- 
scription of finite nuclei and nuclear matter. The present work aims at determining such effective 
interactions for the DDRHF. As already done with the RMF theory j^], the multi-parameter fit- 
ting is performed with the Levenberg-Marquardt method [toI ] . In the DDRHF, 8 parameters for the 
density-dependence of g a and g w are reduced to 3 with the five constraint conditions mentioned in 
Subsec. II. A . Among the effective mesons which mediate the interaction, the pion meson is the most 
important one. Within the RMF model the cr-meson takes into account only the two-pion exchange 
in a phenomenological way. In the Hartree approximation the one-pion exchange does not contribute 
because of parity conservation. This approach cannot take into account the effective spin-spin interac- 
tion coming from the nucleon density through the anti-symmetrization of the wave function. Because 
of its relatively small mass the pion cannot be approximated by a zero range force with a rather trivial 
exchange term. In nuclear matter and in spin-saturated systems the one-pion exchange plays only a 



minor ro 
systems 



e. However, it has an important influence on single-particle properties in spin-unsaturated 



531 ]. In the DDRHF, the contributions of one-pion exchange can be included naturally. 



With similar considerations as for the case of p-N coupling in DDRMF, an exponential density- 
dependent behavior for the 7r-N coupling is adopted here. As further constraints, we impose that the 
coupling constants g p (Q) and f n (0) reach their experimental values at p=0. Thus, there are in total 
8 free parameters in the DDRHF: the mass m a of the a- meson , the couplings constants g a (po) , 
9u){po) an d the 3 parameters describing their density dependence, and the a p and a n parameters for 
the density dependence of the p-N and 7r-N couplings. This gives the new effective interaction PKOl 
shown in Table E where the masses of the nucleon, u-, p- and 7r-mesons are respectively taken as 
M = 938.9MeV, m u = 783.0MeV, m p = 769.0MeV and m n = 138.0MeV, and the experimental values 
for the coupling constants g p and f n are respectively 2.629 and 1.0. In order to investigate the role 
of one-pion exchange, two other parameterizations (see Table QJ are also developed as PK02 without 
7T-N coupling, where g p (0) is free to be adjusted, and PK03 with free g p (0), where a n is adjusted by 
hand. 

In these parameterizations, the masses of the nuclei 16 0, 40 Ca, 48 Ca, 56 Ni, 68 Ni, 90 Zr, 116 Sn, 132 Sn, 



11 



TABLE I: The effective interactions of the DDRHF: PKOl, PK02 and PK03. The masses of nucleon, w-, p- 
and 7r-mesons are respectively taken (in MeV) as M = 938.9, = 783.0, m p = 769.0 and = 138.0. 





m CT (MeV) 


9<y(po) 


9u{po) 


9 P {0) 


/*(0) 


a p 




Po 


PKOl 
PK02 
PK03 


525.769084 
534.461766 
525.667686 


8.833239 
8.920597 
8.895635 


10.729933 
10.550553 
10.802690 


2.629000 
4.068299 
3.832480 


1.000000 
1.000000 


0.076760 
0.631605 
0.635336 


1.231976 
0.934122 


0.151989 
0.151021 
0.153006 




a a 


K 








K 




d u 


PKOl 
PK02 
PK03 


1.384494 
1.375772 
1.244635 


1.513190 
2.064391 
1.566659 


2.296615 
3.052417 
2.074581 


0.380974 
0.330459 
0.400843 


1.403347 
1.451420 
1.245714 


2.008719 
3.574373 
1.645754 


3.046686 
5.478373 
2.177077 


0.330770 
0.246668 
0.391293 



182 Pb, 194 Pb, 208 Pb and 214 Pb are fitted. The bulk properties of symmetric nuclear matter, i.e., the 
saturation density po, the compression modulus K and the symmetry energy J at saturation point 
are also included as constraints to improve the parameterizations. These inputs have been used to 
minimize the least square error, 

N 



X\a) = £ 



(27) 



i=l 

where yi and Oi are the observables and corresponding weights. The vector a is the ensemble of the 
free parameters. Notice that the radii are not included as observables in the fitting procedure. Only 
the masses of the above nuclei and the bulk properties of nuclear matter are taken into account as 
observables. The reason is that the radii are much less sensitive to the changes in the parameters. For 
the parameter search, the best choice is to find a set of observables which have more or less the same 
sensitivity to the parameters. There would be no special difficulty to include the radii in the list of 
fitted data, but we have not made this choice. It is thus more remarkable that the radii calculated 
with our models have a good agreement with the data, as it will be seen in Table V. 

The properties of nuclear matter such as the saturation point, the compression modulus K and the 
symmetry energy J are just empirically determined quantities which are known with some errors and 
therefore, they can be slightly varied to improve the description of finite nuclei. In fact, it is possible 
to obtain fairly good description of radii by properly choosing the values of bulk properties of nuclear 
matter. The isospin related properties in nuclei can also be well described with an appropriate choice 
of the symmetry energy J. In Table[M]are shown the compression modulus K, the saturation baryonic 
density p&, the symmetry ene rgy J, the scalar (Dirac) mass Mg, the non-relativistic and relativistic 
effective masses and |58j ] at the saturation point of nuclear matter. The results are calculated 



by the DDRHF with PKOL PK02 and PK03, the RMF with PK1, PK1R, PKDD JM, NL3 71 1, 
TW99 Q and DD-ME1 3> and the RHF approach with set e in Ref. Q], HFSI [J ZRL1 57]. 
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It is found that the values of the compression modulus K, saturation density po, symmetry energy J 
and the binding energy per particle E/A provided by the DDRHF approach are all in the currently 
acceptable regions [d], namely K ~ [250MeV, 270MeV] , J ~ [32MeV, 36MeV] and E/A ~ 16MeV. 
Compared with earlier RHF applications j^|, the improvement on the compression modulus is quite 
significant. It is worthwhile to mention that the values of effective mass predicted by the DDRHF 
with PKO series are quite close to the empirical value M* = 0.75 0]. 



TABLE II: The bulk properties of symmetric nuclear matter calculated in DDRHF with PKOl, PK02 and 
PK03. The resultspredicted by the RM1 

n n 

with set e in Ref. |55j, HFSI |56j, ZRL1 [5JJ are also given for comparison. 



PK03. The resultspredicted by the RMF theory with PK1, PKDD, NL3 and DD-ME1, and the RHF approach 







^(MeV) 


Po(fm- 3 ) 


J(MeV) 


E/A(McV) 


M* s (p f )/M 




M£/M 


PKOl 


250.24 


0.1520 


34.371 


-15.996 


0.5900 


0.7459 


0.7272 


PK02 


249.60 


0.1510 


32.492 


-16.027 


0.6025 


0.7636 


0.7447 


PK03 


262.47 


0.1530 


32.987 


-16.041 


0.5862 


0.7416 


0.7229 


PK1 


282.69 


0.1482 


37.641 


-16.268 


0.6055 


0.6811 


0.6642 


PK1R 


283.67 


0.1482 


37.831 


-16.274 


0.6052 


0.6812 


0.6639 


TM1 


281.16 


0.1452 


36.892 


-16.263 


0.6344 


0.7074 


0.6900 


NL3 


271.73 


0.1483 


37.416 


-16.250 


0.5950 


0.6720 


0.6547 


PKDD 


262.18 


0.1500 


36.790 


-16.267 


0.5712 


0.6507 


0.6334 


TW99 


240.27 


0.1530 


32.767 


-16.247 


0.5549 


0.6371 


0.6198 


DD-ME1 


244.76 


0.1520 


33.069 


-16.202 


0.5779 


0.6574 


0.6403 


Ref. [55] 


465.00 


0.1484 


28.000 


-15.750 


0.5600 






HFSI 


250.00 


0.1400 


35.000 


-15.750 


0.6100 






ZRL1 


250.00 


0.1550 


35.000 


-16.390 


0.5800 







For the description of finite nuclei, the comparisons between the DDRHF and the previous RHF 
approaches 0, 0, are shown in Table Mill It is found that the DDRHF provides successful 
quantitative descriptions of binding energies, spin-orbit splittings and charge radii for the nuclei 16 O, 
40 Ca, 48 Ca, 90 Zr and 208 Pb , where the deviations from the data on the binding energies are less 
than 1.5 MeV. Compared with the previous RHF results, the present DDRHF brings substantial 
improvements on these observables. 

We now focus on the comparison of results between the DDRHF with PKOl, PK02 and PK03 
and the RMF with the nondinear PK1 NL3 and the density-dependent PKDD Q, DD-ME1 
65l |. The reference nuclei are those selected to determine these seven parameterizations, i.e., 16 0, 
40 Ca, 48 Ca, 56 Ni, 58 Ni, 68 Ni, 90 Zr, 112 Sn, 116 Sn, 124 Sn, 132 Sn, 182 Pb, 194 Pb, 204 Pb, 208 Pb, 214 Pb, and 
210 Po. Table [TV1 shows the masses (binding energies) of these nuclei, where the numbers in boldface 
denote the selected observable used in the corresponding parameterization. The last row of Table ITV1 
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TABLE III: The properties of the nuclei 16 0, 40 Ca, 48 Ca, 90 Zr and 208 Pb predicted by the DDRHF with PK01, 
PK02, and PK03, and the RHF approaches with set c in Ref. Q, HFSI Q, ZRL1 Q The total binding 
energies E, the proton spin-orbit splitting Als for the lp shell of 16 and Id shell of 40 Ca and 48 Ca nuclei are 
given in MeV. The r.m.s. charge rad ii rv are given in fm. The experimental data are taken from Ref. [j^ for 
the binding energies, and from Rcfs. [zJIzflQ for charge radii, and from Ref. jjgl for spin-orbit splitting. 



Set 


16 Q 

E r c A LS 


40 Ca 
E r c A LS 


48 Ca 
E r c A LS 


90 Zr 

E r c 


208p b 

E r c 


EXP 


-127.63 2.693 6.32 


-342.05 3.478 5.94 


-415.99 3.479 5.01 


-783.89 4.270 


-1636.45 5.504 


PKOl 
PK02 
PK03 


-128.33 2.676 6.36 
-126.94 2.659 6.25 
-128.34 2.687 6.20 


-343.28 3.443 6.63 
-340.93 3.427 6.47 
-343.37 3.448 6.63 


-417.37 3.451 5.59 
-415.55 3.450 6.15 
-416.80 3.459 5.27 


-784.61 4.251 
-782.56 4.247 
-784.89 4.254 


-1636.92 5.505 
-1636.63 5.509 
-1637.42 5.505 


Ref. [55] 
HFSI 
ZRL1 


-99.52 2.73 7.3 
-128.64 2.73 6.4 
-127.68 2.71 6.3 


-280.8 3.47 8.0 
-341.2 3.48 7.05 
-341.2 3.44 7.1 


-349.44 3.47 4.1 
-414.24 3.48 3.27 
-417.12 3.49 2.6 


-673.2 4.26 
-779.4 4.26 
-787.5 4.25 


-1406.08 5.47 
-1622.4 5.52 
-1636.96 5.49 



presents the r.m.s. deviations from the data, defined as 

i 

2 

(28) 

It is seen that the DDRHF reproduces well the binding energies of the selected nuclei as compared 
with the data, where few (only 4 cases for PKOl) deviations are larger than 1.5 MeV. From the r.m.s. 
deviations in Table 1 1 VI it is also demonstrated that the DDRHF achieves successful quantitative 
descriptions comparable to the RMF, even slightly better for the selected nuclei. Since the selected 
nuclei cover from the light to heavy regions, one may conclude that the DDRHF with PKO series can 
provide appropriate quantitative descriptions of the masses of finite nuclei throughout the periodic 
table. 

In the DDRHF parameterizations, the radii are not included as fitted observables. The charge radii 
of the selected nuclei calculated by the DDRHF and RMF are shown in Table As one can see from 
the r.m.s. deviations, the DDRHF gives satisfactory quantitative descriptions of the charge radii, less 
good but still comparable to the RMF predictions. It is also found that large deviations only exist in 
the light region. For heavy nuclei, the DDRHF has the same quality of fit as the RMF. It should also 
be kept in mind that other corrections beyond the mean field can affect the radii. For instance, the 
corrections from the RPA correlations may increase the charge radii for 40 Ca. 



A 



N 



E 



EXP 



E 



Cal Y /N 
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TABLE IV: The binding energies of the selected nuclei used in the various models. The results are calculated 
by the DDRHF with PKOl, PK02, and PK03, and the RMF with PKl, PKDD, NL3 and DD-ME1. The 
experimental data are taken from Ref. [7^ . The boldface values correspond to nuclei used in the various fits. 





EXP 


PKOl 


PK02 


PK03 


PKl 


NL3 


PKDD 


DD-ME1 


16Q 


-127.619 


-128.325 


-126.942 


-128.345 


-128.094 


-127.127 


-127.808 


-127.926 


40 Ca 


-342.052 


-343.275 


-340.930 


-343.374 


-342.773 


-341.709 


-342.579 


-343.653 


48 Ca 


-415.991 


-417.372 


-415.554 


-416.799 


-416.077 


-415.377 


-415.944 


-415.012 


56 Ni 


-483.998 


-483.061 


-484.548 


-480.807 


-483.956 


-483.599 


-484.479 


-480.869 


58 Ni 


-506.454 


-502.679 


-503.739 


-500.912 


-504.033 


-503.395 


-504.013 


-501.312 


68 Ni 


-590.430 


-591.370 


-588.581 


-590.693 


-591.685 


-591.456 


-591.241 


-592.253 


90 Zr 


-783.893 


-784.605 


-782.444 


-784.891 


-784.781 


-783.859 


-784.879 


-784.206 


112 Sn 


-953.529 


-951.651 


-950.488 


-951.472 


-954.210 


-952.562 


-953.730 


-952.468 


116 Sn 


-988.681 


-987.276 


-985.808 


-987.119 


-988.491 


-987.699 


-988.066 


-988.470 


124 Sn 


-1049.963 


-1049.619 


-1047.818 


-1048.918 


-1049.162 


-1049.884 


-1048.113 


-1049.880 


132 Sn 


-1102.920 


-1103.585 


-1103.171 


-1102.106 


-1103.503 


-1105.459 


-1102.648 


-1103.857 


182p b 


-1411.650 


-1412.159 


-1411.832 


-1413.027 


-1416.431 


-1415.184 


-1416.309 


-1415.216 


194p b 


-1525.930 


-1523.023 


-1524.502 


-1523.546 


-1525.536 


-1525.733 


-1525.474 


-1524.937 


204p b 


-1607.520 


-1606.269 


-1607.230 


-1606.737 


-1607.851 


-1609.906 


-1607.770 


-1609.676 


208p b 


-1636.446 


-1636.913 


-1636.625 


-1637.435 


-1637.443 


-1640.584 


-1637.387 


-1641.415 


214p b 


-1663.298 


-1660.797 


-1658.841 


-1661.099 


-1659.382 


-1662.551 


-1656.084 


-1662.011 


210p o 


-1645.228 


-1646.350 


-1645.877 


-1647.242 


-1648.443 


-1650.755 


-1648.039 


-1651.482 


A 




1.6177 


1.8745 


2.0489 


1.8825 


2.2506 


2.3620 


2.7561 



TABLE V: Same as Table QV] for charge radii. The experimental data are taken from Refs. 





EXP 


PK01 


PK02 


PK03 


PK1 


NL3 


PKDD 


DD-ME1 

V J \ J 1 V 1 11 1 


16q 


9 693 


9 6763 

Zi.U 1 UO 


9 6^93 

Zi .UOjO 


2 6871 


9 69^7 


2.7251 


9 6988 


9 7968 

Zi . 1 ZiUO 


V_y CI 


3 478 

0.4:10 


3 4498 

O .44ZjO 


3 49fifi 

0.4ZiUU 


3 447Q 

O .44 i c/ 


3 44^ 

.4400 


9 4R7Q 
o .4:U j y 


3 441 & 

0.4:110 


% 4fi99 


48 Ca 


3 4-7Q 


3 4506 

O .40Uvl 


o.^touu 


3 4^04 


3 4^7^ 

O .4U i O 


3 4846 

O .404U 


3 471 6 

O .4 j 1U 


3 4946 

O .4C/4U 


56 Ni 




3 6903 


3 6943 

O . UC/40 


3 6944 

O .Uv/44 


3 7085 


3.7122 


3 71 69 

. 1 lUii 


3 731 5 

O. I O 1. o 


58 Ni 


3 77fi 

O.I IU 


3 71 QQ 

O . 1 ±iJiJ 


3 721 3 

O. 1 1.0 


3 7964 


3 7383 

O . f ooo 


3 7435 

. 1 400 


3 7449 

O . f 44 Zj 


3 761 3 

. 1 UIO 


68 Ni 




3.8501 


3.8518 


3.8613 


3.8621 


3.8773 


3.8681 


3 8926 


90 Zr 


4.270 


4.2505 


4.2466 


4.2537 


4.2522 


4.2689 


4.2534 


4.2725 


112 Sn 


4.593 


4.5635 


4.5640 


4.5672 


4.5704 


4.5861 


4.5722 


4.5901 


116 Sn 


4 69 C I 
4. UZjO 


4 5930 

4. Oi/OU 


4.5944 


4 ^981 

4. Oc/O -1 


4 ^984 

4.0i?04 


4.6149 


4 6004 

4 . UV7V74 


4.6212 


124c 

Oil 


4 677 

4.UI 1 


4.6481 


4 6523 

4. UOZiO 


4 6545 

4.U040 


4 6536 

4.UOOU 


4 6685 

4. UUOO 


4 6567 

4.UOU i 


4 6781 

4.UI O -1 


132 Sn 

Oil 




4 6967 


4 70S3 

4. i UOO 


4 701 1 

4. i \J X ± 


4 7064 

4. 1 UU1 


4 71 83 

4. f lOO 


4 71 09 

4. i ±\J Li 


4 7970 

4. 1 Lt I \J 


182pu 
1 u 




5 371 

O .O ( l-VJ 


5 3687 

O .OUO I 


5 3683 


5 3708 

O.O 1 UO 


5 3900 

O .OJUU 


5 3709 

O.O 1 VI 


5 3896 


i94 Pb 


5.442 


5.4333 


5.4311 


5.4320 


5.4327 


5.4506 


5.4329 


5.4539 


204p b 


5.482 


5.4856 


5.4883 


5.4839 


5.4869 


5.5027 


5.4877 


5.5038 


208p b 


5.504 


5.5054 


5.5093 


5.5046 


5.5048 


5.5204 


5.5053 


5.5224 


214 Pb 


5.559 


5.5633 


5.5644 


5.5629 


5.5658 


5.5820 


5.5635 


5.5779 


210p o 




5.5394 


5.5425 


5.5383 


5.5370 


5.5539 


5.5371 


5.5544 


A 




0.0269 


0.0299 


0.0225 


0.0204 


0.0177 


0.0188 


0.0163 
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IV. CONCLUSIONS 

In this work, we have introduced an effective Lagrangian suitable for a description of nuclear systems 
in a relativistic Hartree-Fock framework. This Lagrangian contains the 4 effective mesons cr,uj,p and 
7r with density-dependent meson-nucleon couplings. The number of adjustable free parameters is 8, 
i.e., comparable with the currently existing Lagrangians where the Fock terms are not treated. We 
recall that the pion can contribute only through its Fock term if parity is to be conserved. The criteria 
for determining the model parameters are the reproduction of the binding energies in a number of 
selected nuclei, and the bulk properties of nuclear matter (saturation point, compression modulus, 
symmetry energy). 

In comparison with other RMF models, the present density-dependent RHF approach is at the same 
level of quantitative description, or slightly better as far as the set of selected nuclei are concerned. This 
applies not only for binding energies but also for observables not considered in the fitting procedure 
such as charge radii. Because the nuclei selected for the parameter fitting cover a wide range of 
masses, it can be expected that the present density-dependent RHF approach can perform quite well 
throughout the periodic table and out to the borderlines of the region of bound nuclei. 

One may wonder if one has gained any benefit over the RMF approach by including the Fock terms. 
The answer is twofold. Firstly, there is less arbitrariness in RHF since one does not drop the exchange 
terms while they are known to be important. Even though their effects can be partly simulated by 
renormalizing the coupling constants, this cannot be completely correct: for instance, the pion will 
never appear explicitly in RMF. Secondly, there are other physical properties not discussed in this 
paper, where the Fock terms and particularly the pion contributions are quite important for explaining 
the nuclear data. This is in particular the case of the isospin dependence of the nucleon effective mass 
[ssT ] , or the evolution of the single-particle energies along isotopic and isotonic chains [3] ■ 

Several aspects remain to be investigated. The next step is to improve the treatment of pairing 
correlations within the RHF by setting up a computational scheme for a relativistic Hartree-Fock- 
Bogoliubov approach. Another extension is to include the tensor coupling of the ^-meson and to 
redetermine the optimal effective Lagrangians for describing finite nuclei. 
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APPENDIX A: THE MICROSCOPIC CENTER-OF-MASS CORRECTIONS 

1. Corrections to the total energy 

The corrections from the center-of-mass motion to the total energy have been systematically dis- 
cussed in Ref. Here for completeness, we just recall the corresponding formula in this subsection. 
The correction of the center-of-mass motion on the binding energy is taken into account in a non- 
relativistic form as, 

where the center-of-mass momentum P c .m. = Yli Pi ano ^ the expectation value of its square is 

(Pirn. ) = Y V %P™ ~ Y V a V bPab ' P*ab + Y V ^ V b u bPab ' Pal • ( A2 ) 
a a,b a,b 

To calculate the above expression, we express P\ m in second quantized form, 



P c.m. = Y P °> ' P b = Y P ^' C a c a' + Y ' P/3/3'44 C /3' C °' ' ( A3 ) 

ab aa' af3a'/3' 



In the spherical case, we obtain 

I 2 J 
.-Pern.) =-Y V ™ Y ( aKm \ A 



anm) 

oik m=—j 

3 J 



Y v <xk v PK {v aK Vf3K + U aK U{3 K ) Y Y \( aKm \ V \P KM )\' 
a,P,K,K rn=—j M=—J 



where 



m=-j r?=± 



and 



2 

J j 1 ' 



flg = J drr^ aK ^ K (A7b) 

In the above expressions, rj denotes the corresponding quantities of the upper (77 = +) or lower (77 = — ) 
components of the Dirac spinor (|18|). 



(A4) 



^ (a«;m| A|cmm) = j 2 Y / drr 2 ^Atp^ , (A5) 



mtl.Mtlj £± [/^LWiJ (A6) 

The corresponding formulas for A, £> read as, 
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2. Corrections to radii 



The definition of the center-of-mass coordinate is 

A 

1 ■ 

r 



~^J2 r i- ( A8 ) 



A 

i=i 



The nucleon r.m.s. radius is defined as (taking protons as representative), 

1 Z 



z 

r 2 

p Z 

i=l 



With the center-of-mass correction, it then becomes 

z 



rl m ., P = \ (01 j>i ~ r cm.) 2 10) = r 2 p + Sr 2 p , (A10) 



i=l 



where the correction term 5r 2 reads as 



2 Z 

6r V = -?E W r * ■ r — 1^) + (01 r c.m. 10) • (AH) 

t=l 



In right hand of above equation, there exist the following type operators 

Similar as the quantization of the center-of-mass momentum, it can be quantized as 

= J ^2(r 2 ) aa >c < a c a > + ^ r aa > ■ rp^c^c^c^ , (A13) 

act 1 aa'/3/3' 

With the BCS state, 

|0) = |BCS) = J] (u k + Ufc 44) |0) , (A14) 

k 

the expectation can be obtained as 

(01 r(i,j) |0) = E^(r 2 ) aQ -E«a^r Q/3 • rp a + ^v a vpu a upr ol p • r^. (A15) 

a a/3 a/3 

Among the right hand of above equation, the terms r a p ■ rg Q and r Q/ g • r 5 s can be obtianed as 

r °> ' r " =S rMm <*>°> > + *« F > F * Ke °" G < J + ® M < • (A16a) 

r a p ■r e $ = j rdrr'dr r £j 2 \^ afi G a Gfi + ^ a pF a F p ] r {d a pG a Gp + ^ a pF a F p ] rl , (A16b) 



where 



* =(-l)H <! f f ! CfcSLO = (-1)^6 J ; J ; ! | Cf° , (Al7a) 



la h 2 ) \ la h 2 
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In the end, the correction terms QA11J1 for the proton and neutron radii can be expressed as, 
H = ~ J [(2rJ " rli) + (»? - r 2 M )] , =~J ~ *1) + - , (A18 



where 



1 f 4 (n 



dr, 



f 2 M =-(Zf 2 p + Nf 2 n ), 



1 Z 

a/3 
1 N 

-TJ ^2 VaV P ( V a v f3 r aP ■ r{3 a - U a Uf3r af3 ■ T^) , 



a/3 



M 



+ -{Zfl + Nf 2 n ) 



a 



(Al9a) 

(Al9b) 
(Al9c) 



3. Corrections to the charge density 



In general, the charge density is obtained from the calculated proton distribution corrected by the 
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. The first correction is done by using the proton 



center-of-mass motion and finite proton size 
density p c .m. in the center-of-mass system, which is related to the HF density through 

4 



PHF 



-^r [e- r ' 2/B2 p,n.(\r-r'\)dri 

B 3 1T2 J 



(A20) 



with B = [h/MivA] 2 . This formula is valid for a system with the harmonic-oscillator center-of-mass 
correction, which gives the center-of-mass correction on the energy as 



E, 



cm. — 



1 



(P 2 ) 
\* cm./ 



2MA 

In term of the center-of-mass momentum, the factor B can be written as, 



b~ 2 = tr l MAuj = h- 2 MAhw = -r 2 {p c 2 ffi ) 

3 



(A21) 



(A22) 



The HF density ()A20|) is then expressed as 



8/2 (PL,) 3 



PHF 



3vr 



exp 



(P 2 V 

\ cm./ 



,'2 



Pc.m. (I r - r'\)dr' . 



(A23) 



Such that, the contributions from the center-of-mass motion to the density can be treated in micro- 
scopic way by introducing the expectation of the center-of-mass momentum (jA3|) . 

The finite proton is also taken into account by convoluting p c .m. with a Gaussian representing the 
proton form factor, 



Pch(r) 



k sm(kr)pnF(k) exp 



\k 2 (B 2 -a 2 ) 



dk , 



(A24) 
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where the Gaussian is (a-^/vr) -3 exp(— r 2 /a 2 ), with a = \/2/3 (r p ) r m , (?>) rms = 0.8 fm and phf(&) 
is the Fourier transform of the HF proton density. When a > B, this expression can be reduced into 
the following form, 

Pch(0 = J dr 2 (a /tt) ~ 3 exp [_( r _ r2 ) 2 / a 2 ] Pc ra ( r ) 

= j dndr 2 (aV^)~ 3 (B^)~ 3 exp [-(r - r 2 ) 2 /a 2 + (r 2 - rrf/B 2 ] pbf(t\) (A25) 

= | dn [(a 2 - 5 2 ) vr] ~ 3/2 exp [-(r - n) 2 / (a 2 - B 2 )] p HF (ri) . 

Denoting by A 2 = l/(a 2 — B 2 ), the integral can be expressed as, 

A 3 



Pch(r) 



TT 



3/2 



J exp [-A 2 (r - r'f] p HF (r')dr' . 



In the spherical case, the integral reads as 
A 3 



Pch(r) 



"vr 3 / 2 
A 3 

: 7T 3 / 2 

A 3 

^3/2 

A 3 



J exp [-A 2 (r - r') 2 ] p U F(r')dr' 

( r' 2 dr'p UF (r') f sintfdtfexp [-A 2 (r - r') 2 ] f 
J Jo Jo 



dip 



j2 dr'p U F(r')exp [-A 2 (r 2 + r' 2 )] / sini?(ii?exp(2A 2 rr'cosi9) x 2vr 



TT 



3/2 



• ,2 drV HF (r')exp [-A 2 (r 2 + r' 2 )] x 



2A 2 



rr 



(A26) 



(A27) 



Finally, one can write the charge density distribution combined with the center-of-mass motion and 
proton size corrections as, 

A 



r'dr'p H F{r') 



e -A 2 (r-r') 2 _ e -A 2 (r+r') 2 



(A28) 
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